AN ASYMPTOTIC LIMIT OF A NAVIER-STOKES SYSTEM WITH 

CAPILLARY EFFECTS 



ANSGAR JUNGEL, CHI-KUN LIN, AND KUNG-CHIEN WU 

Abstract. A combined incompressible and vanishing capillarity limit in the barotropic 
compressible Navier-Stokes equations for smooth solutions is proved. The equations are 
considered on the two-dimensional torus with well prepared initial data. The momentum 
equation contains a rotational term originating from a Coriolis force, a general Kortewcg- 
type tensor modeling capillary effects, and a density-dependent viscosity. The limiting 
model is the viscous quasi-geostrophic equation for the "rotated" velocity potential. The 
proof of the singular limit is based on the modulated energy method with a careful choice 
of the correction terms. 



I. Introduction 

The aim of this paper is to prove a combined incompressible and vanishing capillarity 
limit for a two-dimensional Navier-Stokes-Korteweg system, leading to the viscous quasi- 
geostrophic equation. We consider the (dimensionless) mass and momentum equations for 
the particle density p(x,t) and the mean velocity u(x, t) = (ui(x,t),U2(x,t)) of a fluid in 
the two-dimensional torus T 2 : 

(1) d t p + div(pw) =0 in T 2 , t > 0, 

(2) d t {pu) + div(pu ® u) + pu L + Vp(p) = div(i\~ + S), 



with initial conditions 



p(.,0) = p°, u(-,0) = u° inT 2 . 
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Here, pu 1 - describes the Coriolis force, u 1 - = (—u 2 , the function p(p) = p 7 /7 with 7 > 1 
denotes the pressure of an ideal gas obeying Boyle's law, K is the Korteweg-type tension 
tensor and 5" the viscous stress tensor. 

More precisely, the free surface tension tensor is given by 

div if = K pV(a'(p)Aa(p)), 

where k > 0, which can be written in conservative form as 

(3) div if = Ko div ((aS(p) - V(p)|V/f) I - Va(p) Va(p)^ , 

where S'(p) = pa'(p) 2 , a(p) is a (nonlinear) function, and I denotes the unit matrix in 
M 2x2 . For a general introduction and the physical background of Navier-Stokes-Korteweg 
systems, we refer to [7J [HJ [19]. In standard Korteweg models, n{p) = u'(p) 2 defines the 
capillarity coefficient [HI Formula (1.29)]. In the shallow-water equation, often a(p) = p 
is used such that div A' = pVAp (see, e.g., [5J[25]). Bresch and Desjardins [5] employed 
general functions a(p) and suitable viscosities allowing for additional energy estimates (also 
see [20]). If a(p) = y/~p, the third-order term can be interpreted as a quantum correction, 
and system ([I])-© (without the rotational term) corresponds to the so-called quantum 
Navier-Stokes model, derived in [Hj and analyzed in [19] . 
The viscous stress tensor is defined by 

div 5 = 2 div (fi(p)D(u)), 

where D(u) = \(Vu + Vu T ) and p(p) denotes the density- dependent viscosity. Often, 
the viscosity in the Navier-Stokes model is assumed to be constant for the mathematical 
analysis [13]. Density-dependent viscosities of the form p(p) = p were chosen in [5J and 
were derived, in the context of the quantum Navier-Stokes model, in [8]. The choice 
/i(p) = cr(p) allows one to exploit a certain entropy structure of the system [6]. 

Without capillary effects, system ([I])-© reduces to the viscous shallow-water or vis- 
cous Saint- Venant equations, whose inviscid version was introduced in [27]. The viscous 
model was formally derived from the three-dimensional Navier-Stokes equations with a 
free moving boundary condition [H]. This derivation was generalized later to varying 
river topologies [25] . The existence of global weak or strong solutions to the Korteweg- 
type shallow-water equations was proved in [61 [7J, [T5J, [T6], [18] under various assumptions 
on the nonlinear functions. In [7J, the authors obtained several existence results of weak 
solutions under various assumptions concerning the density dependency of the coefficients. 
The notion of weak solution involves test functions depending on the density; this allows 
one to circumvent the vacuum problem. Duan et al. [TU] showed the existence of local 
classical solutions to the shallow-water model without capillary effects. For more details 
and references on the shallow- water system, we refer to the review [1]. 

The combined incompressible and vanishing capillarity limit studied in this work is based 
on the scaling t et, u \-> eu, p(p) \-> ep>(p) and on the choice n$ = e 2a (0 < a,e < 1), 
which gives 

(4) 8 tPe + dw{p e u e ) = in T 2 , t > 0, 
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(5) d t {p £ u £ ) + dw(p £ u £ g> u £ ) + -p £ uj + -^V(pJ) - 2e 2 ^p £ V(a'(p £ ) Aa(p £ )) 

= 2div(p(p £ )D(u £ )), 
with the initial conditions 

(6) p £ (-,0) = p°, u £ (;0)=u° £ inT 2 . 

The condition a < 1 is needed to control the capillary energy; see the energy identity in 
Lemma [2] below. The local existence of smooth solutions to (ED)-© is discussed in Appendix 

m 

When letting e — > 0, it holds p £ — > 1 and p £ u £ — > V^cfi = (— <90/<9x 2 , <90/ctei) in 
appropriate function spaces, where solves the viscous quasi-geostrophic equation [221 
Chapter 6] (see Section [2] for details) 

(7) <9 4 (A0 - 4>) + (V^ • V)(A0) = /i(l)A 2 in T 2 , t > 0, 

(8) <f>(-,0) = <f>° inT 2 . 

The objective of this paper is to make this limit rigorous. Our proof requires the (local) 
existence of a smooth solution to ©-([8]), which is shown in Appendix [A] Several derivations 
of inviscid quasi-geostrophic equations have been published; see, e.g., [HJCI2JI2E]- The reader 
is also referred to the monograph [23] for a more complete discussion of this model. The 
viscous equation was derived rigorously for weak solutions from the shallow-water system in 
[5]. The proof is essentially based on the presence of the additional viscous part div(pVn) 
and a friction term in the momentum equation. The novelty of the present paper is that 
these expressions are not needed and that more general expressions can be considered. 
In particular, we allow for viscous terms of the type div (p(p)D(u)), and no friction is 
prescribed. 

In the following, we describe our main result. In order to simplify the presentation, we 
assume that the nonlinearities are given by power-law functions: 

a(p)=p s , p(p) = p m forp>0, 

where s > and m > 0. The exponents s and m cannot be chosen freely; we need to 
suppose that 

1 T + 1 

(9) 0< s < 1, m = s + -<- . 

K J - > 2 ~ 2 

This assumption includes the quantum Navier-Stokes model s = 1/2, m = 1 and the 
shallow- water model with s — 1, m — 3/2. Furthermore, we assume that the initial data 
are sufficiently regular (ensuring the local-in-time existence of smooth solutions) 

p o e H k(j2^ u o e ffk-i^ ^ e ffk+i^ where k > 2) 
and that they are well prepared: 

(10) C7 £ (0°) 0°, e-\pl - 1) 0°, ^ £ u° £ VV, e^Vy/fl -+ 
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in L 2 (T 2 ) as e ->■ 0, where p £ = 1 + e</>° (this defines </> £ ), 

(11) G £ (0 £ ) = ^sign(0 £ )V/i(l + £0 e ), p £ = 1 + e0 £ , 

and the internal energy /i(p) is defined by h"(p) = p\p)j p = p 7 " 2 and h(l) = h'(l) = 
(see (fT3|) for an explicit expression). Note that the convergence £ -1 (p £ — 1) — > <p° in L 2 (T 2 ) 
implies that G £ {<p° £ ) 4>° in L X (T 2 ) if p° is bounded in L°°(T 2 ) (see $17])). 

Theorem 1. Zet < a < 1 and 7 > 1. We suppose that (Q holds and that the initial 
data satisfy ffTUl) . Furthermore, let (p e ,u e ) be the classical solution to (SJ) -(JSJ) and let <fi be 
the classical solution to (JTJ) -(JHD , both on the time interval (0,T). Then, as e — > 0, 

p £ -^l znL°°(0,r;L 7 (T 2 )), 

p e u £ -)• V x L°°(0, T; L 27 /^(T 2 )). 

Furthermore, if s < ~ and 7 > 2(1 — s) or if s = 1 and 7 > 2, 

p £ ^l mL°°(0,T;L p (T 2 )), 

PA^V 1 ^ znL^^TjL^T 2 )), 

/or all 1 < p < 00 and 1 < q < 2. 

The proof is based on the modulated energy method, first introduced by Brenier in a 
kinetic context [2] and later extended to various models, e.g. [U El [22]. The idea of the 
method is to estimate, through its time derivative, a suitable modification of the energy 
by introducing in the energy the solution of the limit equation. We suggest the following 
form of the modulated energy: 

H e (t) = j j2 (fk - VVl 2 + \\G e {<Pe) ~4>\ 2 + 2e 2(a - 1) |W(P £ )| 2 ) dx 

(12) +2/ / p(Ps)\D(u £ ) - D{V^4>)\ 2 dx, 

Jo Jt 2 

These terms express the differences of the kinetic, internal, and Korteweg energies as well 
as the viscosity. Differentiating the modulated energy with respect to time and employing 
the evolution equations, elaborated computations lead to the inequality 

H £ {t)<c[ H E (s)ds + o(l), t>0, 
Jo 

where o(l) denotes terms vanishing in the limit e — > 0, uniformly in time. The Gronwall 
lemma then implies the result. 

The paper is organized as follows. In Section [2J we derive the energy identities for the 
shallow-water system and the quasi-geostrophic equation and give a formal derivation of 
the latter model from the former one. Theorem [1] is proved in Section [3j In the appendix, 
we discuss the existence of local smooth solutions to (J4j)— (J5j) and give an existence proof 
for local smooth solutions to ([7])- ([8]). 
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2. Auxiliary results 

In this section, we derive the energy estimates for (Hj)-© and derive formally the quasi- 
geostrophic equation (J7J). Based on the definition h"(p) = p'(p)/p, h(l) = h'(l) = 0, we 
can give an explicit formula for this function: 

(13) ^(P) = 7(7 1 _ 1) (p 7 -l-7(p~l)), P>0. 

The energy identity for (ED)-© is given as follows. 

Lemma 2. Let (p £ ,u £ ) be a smooth solution to (jl])-© on (0,T). Then the energy identity 

dE F „ , „. 

-£ + D e = 0, te(0,T), 

holds, where the energy E £ and energy dissipation D £ are defined by, respectively, 
E *= f (\h{p £ ) + \p £ \u £ \ 2 + 2e 2 ^\Va{p £ )\ 2 ]dx, D £ = 2 [ p{p £ )\D{u £ )\ 2 dx. 

Proof. Multiply (0} by e~ 2 h'(p e ) - §|u £ | 2 - 2e 2{a ~ 1) a' (p £ ) Aa(p £ ), integrate over T 2 , and 
then integrate by parts: 

= / (^dth(p e ) - \h"(p £ )Vp £ ■ {p £ u £ ) - \ \u £ \ 2 d t p £ + p £ u £ ■ Vu £ ■ u £ 

+ 4s 2 ^Va(p £ ) ■ Vd t a(p E ) - 2e 2 ^ div(p e u e y (p e )A<7(p e ))dz. 
Multiplying ([5]) by u £ and integrating over T 2 gives, since ujr ■ u £ = 0, 

= / (dt(pu £ ) ■ u £ - p £ [u £ ® u £ ) : Vu £ + \p]^Vp £ ■ u £ 
Jx 2 v £ 

+ 2s 2 ^a'(p £ )Aa(p £ ) div(p £ u £ ) - 2p{p e )D{u e ) : Vu £ )dx, 



where ":" means summation over both matrix indices. Observing that h satisfies h"(p £ ) = 
p1~ 2 and using the identity D{u £ ) : V« e = |-D(u e )| 2 , the sum of the above two equations 



becomes 



d f ( 1 , , , 1 



dtU \e^ ) + 2 p ^ U 



2e 2 ^\Va{p £ )\ 2 ]dx + 2 [ p(p £ )\D(u £ )\ 2 dx = 0, 
/ Jt 2 



which proves the lemma. □ 

A consequence of the energy identity is the following estimate. 

Lemma 3. Let (p £ ,u £ ) be a smooth solution to (SJ) -(EJ) on (0,T). Then there exists C > 
such that for all < e < 1, 

(14) \\ P e ~ l|U°o(0,T;L,(T 2 )) < Ce"**^™ if! > 1, 

(15) \\p £ - 1||l-(o,t ; ,l 2 (t 2 )) < Ce if^>2. 
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Proof. If 7 = 2, h(p) = |(p — l) 2 , and the result follows immediately from Lemma [2j Let 

7 > 2. We claim that h(p) > \p — 1| 7 / (7(7 — 1)) for p > 0. Then the result follows again 
from the energy identity. Indeed, the function f(p) = p 7 — 1 — j(p — 1) — \p — 1| 7 is convex 
in (|, 00) and concave in (0, ~). Since the values /(0) = 7 — 2 and /(|) = 7/2 — 1 are 
positive, / > on [0, |]. Furthermore, /(l) = /'(l) = which implies, together with the 
convexity, that / > in [|,oo), proving the claim. Finally, let 7 < 2. By [231 P- 591], 
h{p) > c R \p — 1| 2 for p < R and h(p) > cr\p — 1| 7 for p > R, for some cr > and R > 0. 
Hence, using Holder's inequality and 7 < 2, 

Pe ~~ 1|Il7(t2) ^ C 

<C[ f h(p £ )dx] +C f h(p £ )dx < C(e 7 + e 2 ) < Ce\ 

where here and in the following C > denotes a generic constant not depending on e. 
Estimate (I15p for 7 > 2 follows from 



X 7/2 , 




\ Pe - l\ 2 dx + / 


\Pe~l\ 



lp £ -l" 2 



L 2 (T 2 ) 



/ (p e - l) 2 dx <C [ h(p £ )dx < Ce 2 

Jj 2 Jj 2 



which finishes the proof. □ 

We perform the formal limit e — > in (Pil)-([5]). For this, we observe that (jlj) can be 
written in terms of <fi £ = (p £ — l)/s as follows: 

d t <f) £ + div(0 e M e ) + - div w £ = 0. 

We apply the operator div -1 (defined by div" L (f 1 ,f 2 ) = —dvi/dx 2 + dv 2 /dxi) to (jSJ) and 
observe that div -1 (p £ ujr)/s = divu £ /e + div(0 e w e ) = —d t (j) £ , by the above equation. Then 
we find that 

d t div ± (p e w e ) + div 1 - div(p £ u £ <g> u £ ) - d t <j> e 

(16) = 2e 2 ^div ± (p £ V(a'(p £ )Aa(p £ ))) + 2div ± dw(fi(p £ )D(u £ )). 

By the energy estimate, p £ — > 1 (in L°°(0, T; L 7 (T 2 ))). Assuming that (f> £ — > (f> and u £ — > 
V" 1 "^ in suitable function spaces and employing the relations 

div x div(V x <8> V x 0) = (V x • V)(A0), 2 div x div( J D(V ± 0)) = A 2 0, 

the formal limit in ( |T6l) yields the limit equation ((Zj). The initial condition reads as 0(-, 0) = 
0°, where 0° = lim e _>o e (-, 0) in T 2 . The energy and the energy dissipation of (ED) equal 

E = I [ (|V0| 2 + 2 )cfe, D = 2^(1) / \D{V^(l>)\ 2 dx. 



2 _ 

Multiplying the limiting equation by and using the properties 

/ (V x 0- V)(A0)0cfe = 0, / (A0) 2 cfe = 2 / ^(V^)! 2 ^, 
Jt 2 Jt 2 Jt 2 
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we find the energy identity of the viscous quasi-geostrophic equation: 

dE 



dt 



D = 0, t>0. 



3. Proof of Theorem Q] 
First, we prove the following lemma. 
Lemma 4. Let T > 0, 7 > 1, and < a < 1. Then 

lim H £ (t) = uniformly in (0,T), 

where H £ is defined in (1121) . 

Proof. Using the definitions of the energy and energy dissipation as well as the relation 

e~ 2 h{p £ 



\G E {<j> e f 



we write 

-t 



H £ (t) = (E £ + E)(t) + / (D £ + D)(s)ds + I f (p £ -l)|V J 

Jo 1 Jt 2 



2 dx 



T 2 

+ 2 / / (jjl{p, 
Jo Jt 2 

h + • • • + h. 



4> £ )4>dx — j p £ u £ ■ V cj)dx - 
Jt 2 

) - i2(l))\D(V L (j))\ 2 dxds - 4 



cj) £ (j)dx 



Jt 2 



p(p £ )D(u £ ) : D{V ± (j))dxds 



The aim is to estimate dH £ /dt. To this end, we treat the integrals Ij or their derivatives 
term by term. By the energy estimates, -^(h + I2) = 0. The integral -Z3 cancels with a 
contribution originating from J 5 ; see below. The estimate of I4, . . . , Jg (or their derivatives) 
is performed in several steps. 

Step 1: estimate of I4. L'Hopital's rule shows that for 7 > 1, 

>? 2 ' ^0 z V z 2 2 



lim 
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Therefore, there exists a nonnegative function /, defined on [0, 00), such that h(l + z) — 
^z 2 f(z) for z > 0, and a function g, defined on [0, 00), such that f(z) — 1 = zg{z) for z > 0. 

Furthermore, the inequalities f(z) > f(0) = 1 and 1^(^)1 < C(l 



,(7-3)^ 



hold, where 



max{0,z}. Finally, we claim that f(z) = 2h(l + z)/z 2 > 2(1 + z) T_2 /(7(7 - 1)) for 
z > and 7 > 4. Indeed, the function w(z) = h(l + z) — z 2 (l + z) 7-2 / (7(7 — 1)) is convex 
in [0, 00) and w(0) = w'(0) = 0, which implies that w(z) > in [0, 00), proving the claim. 
With these preparations, we can estimate the difference G £ ((f> £ ) — <p £ appearing in J 4 : 



1 a 




10,1 ^^J{^7)-l 
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In view of the bounds for / and g as well as the relation e<p £ — p e — 1, we infer that 

r 1 + n (7_3) 

(17) |G £ (0 e )-0 £ |<-| Pi 



+ 



This bound allows us to estimate I4. Indeed, if 1 < 7 < 4, by (fl4j) . 

< jUh^L^WPe - l||ic (0iT;Ll(T2)) < Ce 2 ^ 1 ^" 1 = o(l) 

uniformly in (0, T). Here and in the following, the constant C > depends on and its 
derivatives but not on e. If 7 > 4, we have, using the upper bound of f(z) for 7 > 4, ( TIT]) . 

and 1 + E(fi £ = p £ , 

\Gs(<i>s) - m < % - - < f b £ - ii 2 (i 

£ Cp e + 1 £ 

We employ estimates (TT4]) -( p~5|) and Holder's inequality to conclude that 

hit) < C||0|| ioO ( O)T . L oo( T 2))£ _1 ||p e — l||L°°(0,T;L 2 (T 2 ))||Pe — 1|U°°(0,T;LT(T 2 )) 
X (l + IIPellL°°(0,T;L^(T 2 ))) 

< Ce^ll^Hioo^T;^^ 2 )) = o(l). 

S'tep 2: estimate of dl^/dt. Inserting the momentum equation (jSJ) and integrating by 
parts, it follows that 

-jr = - [ d t (p £ u £ ) ■ V^^dx - I p £ u £ ■ V L d t 4>dx 
at Jj2 ./ ••' 

p £ (u £ <S> u £ ) : W ± (f>dx + - / p £ uf ■ V~ L (j)dx 

£ Jt 2 

+ -5-/ / p e ^(a , {p e )Aa(j) e ))-V- 

£ 7 Jt 2 </t 2 

+ 2 / n(p e )D(u e ) : W^cfe - / p e it e ■ V^d^dx 
= Ji + • • • + J 6 . 

We treat the integrals Ji, . . . , Jq term by term. The integral J2 can be written as 

■h = - I p £ u £ ■ V<f)dx. 
e Jt 2 

The third integral vanishes since div V ± = 0: 

J 3 = -^ / Pi div(V ± 0)dx = 0. 
£ 2 7 Jt 2 

Using the identity fl3]) and div V -1 = 0, we compute 



J A = e 2 ^ 



(AS(p e )- V(p £ )|Vp,| 2 ) div(V J 
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- (Vcr(p £ ) <g> V«7(p B )) : VV 1 ^! 
<C# £ . 

Integration by parts and using div V 1 " = again yields 
J 5 = - / M/°eK • (V^A^ + V div(V ± <p))dx 

JT 2 

-I p'(p E )(Vp E ®u E + u £ ®Vp E ) :VV ± <pdx 
Jt 2 

= - n{Pe)u e ■ V ± A4>dx 

JT 2 

- 2 f £} Pe / J Va(p £ ) ® {^p~u £ ) + (v^n £ ) <g> Va(p £ )) : VV^ds. 

The assumptions on /x and a (see (jHJ)) yield p'(p £ )/ (y/p^a'(p e )) = pT s ■ Hence, applying 
the Cauchy-Schwarz inequality, the last integral is bounded from above by 

CUVa^H^^ll^ell^^) < Ce 2 ^ = o(l). 

We conclude that 

J 5 <~ MPeK- V ± A0rfx + o(l). 
Jt 2 

The integral Jq remains unchanged. Finally, we estimate J\. To this end, we add and 
substract the expression V" 1 ^ such that J\ — K\ H — ■ + K4, where 



Ki = - / pe{u E - V^0) <g> [u e - V^0) : VV^cfe, 



K 2 = ~ / p £ V x ® u £ : VV^cb, 

JT 2 

K 3 = - / p £ w £ (8) V ± <p : VV^dx, 



AT 4 = / p £ V^0 ® : VV^(f)dx. 
Jt 2 

The first integral can be bounded by the modulated energy: 

K\ < C I p E \u E - V ± 4>\ 2 dx < CH E . 
Jt 2 

A reformulation yields 

K 2 = - [ p E u E ■ {(y ± <p ■ W)W ± <P)dx. 

JT 2 

We employ the continuity equation (@J to find 

K 3 = -- / Pe u £ ■ V\V ± 4>\ 2 dxd = - I div(p £ u £ )\V ± 0\ 2 dx 

I Jj2 I Jj2 



10 A. JUNGEL, C.-K. LIN, AND K.-C. WU 



\ [ d t (p £ -l)\V L <P\ 2 dx 
1 Jt 2 

[ (p £ -l)\V^\ 2 dx + \ [ (p £ -l)d t \V L <P\ 2 dx 
Jt 2 z Jt 2 



2 

1 d 

2dt 



ix 



Finally, using again div V -1 = 0, 

K A = - [ p e ((V ± • V)V ± 0) • V^dx 
Jt 2 

= - [ (pe-1) ((VV • V)V ± 0) • V^dx - [ ((V x • V) V x 0) • V- 

Jt 2 

= - [ (Pe-l) ((V^ • V) V ± 0) ■ V X (j>dx - \ I V x • V(| V L <P\ 2 )dx 
Jt 2 * Jt 2 

= - [ (p £ -l)((V ± 0-V)V ± 0) ■ V x 0rfx + ^ f div(V ± 0)|V ± < />r^ 
J T 2 2 Jt 2 

= o(l). 

In the last step, we have employed estimate (fT4l) for p e — l. Summarizing the estimates for 
Ki, . . . , K4, we have shown that 

■h < CH e - ^ - / ((V ± • V)V x 0) • {p £ u £ )dx + o(l). 
at Jj2 

Then, summarizing the estimates for J 1; . . . , J 6 , we obtain 

^ < CH £ - ^ + - [ p £ u £ ■ Vc^dx 
dt dt e Jf 2 

-I {(d t + V x ■ V) V x + /i(l) V x A0) • (p £ u £ )dx 
Jt 2 

- I (p(p £ )-p(l)p £ )u £ - V ± A0c/x + o(l). 

The last integral can be estimated by employing the assumptions on p and Holder's in- 
equality: 

/ AHPe) ~~ p{^)Pe , — n l A , , ^ n\\ to- 1/2 1/211 || / — II 

/ p= yfp £ u £ -V A<pdx<C\\p £ i - pj ||l 2 (t 2 )|Ivp7w £ ||l 2 (t 2 )- 

7T 2 VPs 

We claim that the first factor on the right-hand side is of order o(l). To prove this 
statement, we consider first ~ < m < 1: 

11 m-l/2 1/2M2 < / n 2m-li _ 1 |2(l-m)j < 11 ||2r7i-l || _1|| 2 (1-™) 

HPe Pe IIl 2 (T 2 ) ^ / Pe \Pe 1 \ aX IIPe|ll,7(T 2 ) HPe 1 1 1 Z,p(T 2 ) ' 



where p = 2j(l — m) / (j — 2m +1). The inequality p < 7 is equivalent to 7 > 1. Note that 
the Holder inequality can be applied since we supposed that 2m — 1 < 7; see ©. Second, 
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let 1 < m < 2 (the case m — 1 being trivial). We compute 

II m-l/2 1/2 1 1 2 ^ I | i |2(m-l) j ^ || || || i ||2(m-l) 

Ik ' ~PJ Lapp) < / Pe|p £ -l| l < IIP £ ||L7(T2)||p e - 1|| L , (T 2) ' 

Vt 2 

where g = 2 / ~f(m — l)/(7 — 1), and g < 7 if and only if m < (7 + l)/2. Finally, if 
2 < m < (7 + l)/2, we find that 

llpr V2 " P £ /2 |l! 2 (T 2 ) <cj Peil + Pr 2 ) 2 \Pe - M 2 dx < C(l + INI^yik - l||i, (T2) , 
with r = 27/(7 — + 3) satisfying r < 7 if and only if m < (7 + l)/2. We conclude that 

for some > 0, and together with ( fl4|) . this shows that the integral is of order o(l). 
Therefore, 

rfJs . ^ w dl 3 1 [ 

ctt at e Jj2 

(18) -/ ((d t + V x 0- V)V ± + /i(l)V ± A0) • (p e u e )cte + o(l). 

JT 2 



S'tep 5 1 .' estimate of dle/dt. Employing (j4j) and ([7]), we can write 

die 
dt 



^ = - / d t (f) £ (f)dx - I (j) £ d t (f)dx 
It 2 Vt 2 



= -/ div(p £ u £ )(j)dx - [ ((d t + V X (f>- V)(A0)-/i(l)A 2 0)0 £ cb 
e Jt 2 </t 2 

(19) =--/ p £ u £ -V(pdx+ I ((d t + V x 0- V)(V x 0) -p(l)V x A0) • V L (j> £ dx. 
£ Jt 2 Jt 2 

We observe that the first integral on the right-hand side cancels with the corresponding 
integral in ({T8l) . To deal with the second integral, we employ again the momentum equation 
©. We write 

-Vp< = (7 - l)Vh(p E ) + V(p £ - 1) = (7 - 1) V% £ ) + eV0 £ . 

7 



Then, because of (u^r) 1 - = —u £ , §5§ is equivalent to 

V x e = p £ u £ - eF £ L , 

where 

7 - 1 

F £ = d t (p £ u £ ) + div(p £ u £ <S> u £ ) H —Vh(p £ ) - 2 div(//(p e )D(u e )) 



£ 2(a-l) 



(vAS^) - ^(^'(pJIp.l 2 ) - div (Va( Pe ) <g> Va(p £ 
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Replacing V -1 ^ in the second integral in f|T9|) by the above expression gives 

/ ({d t + ■ V)(V x 0) - A*(1)V- L A0) • V^dx 
Jt 2 

= [ {(dt + V x • V)(V x 0) - /i(l)V ± A0) • (p £ u e - 
./-••* 

We claim that the integral containing Fj- is bounded in an appropriate space. Indeed, 
let ip be a smooth (vector-valued) test function. The first term of F £ is written in weak 
form as follows: 

/ d t (p £ u £ ) ■ ipdxds = - / p £ u £ ■ d t ipdxds + / (p £ u £ )(t) ■ ijj(t)dx 
o Jt 2 Jo jt 2 Jt 2 



p° £ u° £ ■ ip(0)dx. 



S«2 

T 2 



These integrals are bounded if p £ u £ is bounded in L°°(0,T; L 1 (T 2 )). This is the case, since 
mass conservation and the energy estimate show that 

\p £ u £ \dx < - I p £ dx H — / p £ \u £ \ 2 dx 



is uniformly bounded in (0,T). An integration by parts gives 

div(p e u £ ® u £ ) ■ ipdxds — — / p £ u £ £g> u £ : Vipdxds, 
'o Jt 2 Jo Jj 2 

and this integral is uniformly bounded, by the energy estimate. Furthermore, again inte- 
grating by parts, 

^-^Vh{p £ ) -2div(p(p £ )D(u £ )) ) ■ i)dxds 



o Jt 



s 

r-T 



^^h(p £ )I-2p(p £ )D(u £ ) ) : Vt/jdxds, 
lo Jt 2 V e J 

which is uniformly bounded since we can estimate 

/ / \p,{p e )D{u E )\dxds < - / / n(p e )dxds + - / / p(p £ )\D{u £ )\ 2 dxds 
Jo Jt 2 2 J Jj2 l J Jj2 

and p(p £ ) < C(l + pi). Also the remaining terms are bounded since 



£2(0.-1) 



jf (vA(S(p e ) - 5(1)) - ^V(S"(p £ )|Vp £ | 2 ) - div(V<7(p £ ) <g> W(p £ ))) • ^cbcfe 
r / ((5(p e )-5(l))Adiv^ + V(p £ )|Vp £ | 2 div^ 



r-T 

(Vcr(p e ) (8) Vcr(p e )) : Vip)dxds. 
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Using the Holder continuity of S(z) = (s/2)z 2s , z > 0, the first summand can be estimated 
by C|p £ — l| mm ^ 1,2s ^. We infer that the corresponding integral is of order o(l). We formulate 
the second summand as 

\e^ a - 1 \2s - I) f [ \Va(p e )\ 2 dwtpdxds. 
2 Jo Jj 2 

In view of the energy estimate, this integral as well as the third summand are uniformly 
bounded. This shows that 

/ ((d t + V ± <j> ■ V)(V x 0) - p^V^) • V ± <f) £ dx 
Jt 2 

= / ((9 t + V ± 0-V)(V ± 0)- yU (l)V ± A0)-(p £ n e )dx + o(l), 
Jt 2 



and consequently, (1191) becomes 
dh 1 



p £ u £ ■ Vcf)dx 

dt s 

+ / ((dt + V x 0- V)(V ± 0)-M1)V ± A0) .( PeUe )dx + o(l). 

Step 4-' estimate of dl-j/dt. The function p satisfies \p(z) — p(l)| = \z m — 1| < \z — l| m 
if m < 1 and \p(z) - p(l)| < C(l + z rn ~ 1 )\z - 1| if m > 1, for z > 0. Therefore, if m < 1, 
taking into account (1141) . 

^ 7 <Ol| n _ 1 ||m II n/'V7 ± ^'l II 2 / mmin{l,2/7} 

j, ^ *\\Pe 1 llL°°(0,T;LT(T 2 )) ll-^l V V) llL°°(0,T;L 2 ^/(^- m )(T 2 )) — ° fc 

Moreover, if 1 < m < (7 + l)/2, using Holder's inequality, 

—j^- < C(l + ^PeIL^^T;^™- 1 )^-^^ 2 ))) 1 1 Per — 1 ||l°°(0,T;L7(T 2 )) < (75 mln { 1 > 2 /l'} . 

The norm of p £ is uniformly bounded since (m — 1)7/ (7 — 1) < 7 is equivalent to m < 7. 
Step 5: estimate of dl 8 /dt. Integration by parts yields 

[ p'(p e ) Vp e ® m £ : VV-^dx + 2 / p(p £ )w e • V ± A(f)dx 
dt ./ Jt 2 

T 2 ^p~a\p £ ) 

+ 2 (p(p £ ) - p(l)p £ )u £ ■ V ± A<f)dx + 2p(l) / p £ u £ ■ V ± A<j)dx. 
By definition of p and a (see ©), it follows that 

^ < C'||Va(p e )|| L 2 (T 2 ) || v /^U e || L 2 (T 2 ) + C||p^~ 1/2 - pl /2 \\L2(T2)\\T/plUs\\L2(T2) 

+ 2p(l) / p £ u £ ■ V ± A(j)dx. 
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Because of the energy estimate, the first summand is of order o(l). The second summand 
has been estimated in Step 2, and it has been found that it is also of order o(l). This 
shows that 



-l<2/i(l)/ p £ u £ - V ± A^dx + o(l). 
at Jj2 



Step 6: conclusion. Adding the estimates for dl^/dt, . . . ,dl 8 /dt, most of the integrals 
cancel, and we end up with 

Integrating over (0, t) gives 

H £ (t) < H e (0) + C [ H £ (s)ds + hit) - J 4 (0) + o(l). 
Jo 

By Step 1, hit) = o(l). Furthermore, h(0) = o(l) by assumption. It holds that H e (Q) = 
o(l) since 

II - vV)IIl»cp) ^ llvW - v ± 0°|U 2( T 2 ) + ||(i - VSvVIIlw 

< || x/pf^ - v^/n^^) + ||i - p 0|| i2(T2) ||v x ^°|U-(T2) 

= 0(1) 

and since the initial data are well prepared. Then the Gronwall lemma implies that H e (t) = 
o(l) finishing the proof. □ 

We are now in the position to prove Theorem [T] which is a consequence of Lemma HI 
We observe that by (JHJ), p e — > 1 in L°°(0,T; L 7 (T 2 )) and, using the Holder inequality and 

27/(7 + 1) < 7, 

\\p £ U £ - V" L 0|| LO o( OiT;i 2 7 /( 7+ i)( T 2) ) < ||v / Pr|U° o (0,r;L27(T 2 ))||\/Pe( U e _ ^ ±< l ) ) II L°°(0,T;L 2 (T 2 )) 

+ I l^e ~~ 1||l°°(O,T;L 2 ^/^+ 1 )(T 2 ))I|V" L 0||l«>(O,T;L°°(T 2 )) 

(20) < C\\yTp- £ (u £ - V- L 0)|| £ ao (OiT . i 2( 1 «)) 

+ C||Pe — 1|U°°(0,T;LT(T 2 ))- 

We conclude that p £ u e -»■ in L°°(0,T; L 2 ^ +1 )(T 2 )). 

Next, let 7 > 2(1 — s) and < s < 1/2. Because of assumption (jHJ), i.e. 7 > 2s, we have 
27/ (7 + 2(1 — s)) < 7, and hence, 

p £ -> 1 in L°°(0,T;L 27/(7+2(1 - s)) (T 2 )) 

as £ — j- 0. Furthermore, since a < 1, Vcr(p £ ) — )■ in L°°(0,T; L 2 (T 2 )) as e — > and thus, 
by Holder's inequality, 

(21) ||V(p £ -l)|| 

L oo (0,T;L 2 ^/(7+ 2 ( 1 - s ))(T 2 )) ~~ L oo (0,T;L 2 7/h+ 2 ( 1 - 3 ))(T 2 )) 

- IIPellL^(0,T;L7(T 2 ))ll V(T (^)lli oo (0,r;L 2 (T 2 )) 0. 
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We infer that p £ -> 1 in L°°(0,T; l^WMH)(f)), Because of the continuous embed- 
ding iyi,27/(7+2(i- s ))( T 2) ^ £7/(i-«) (T 2 ), this implies that p £ -> 1 in L°°(0,T; L 7 /( 1_S )(T 2 )). 
Since 27/(7 + 2 (! - s) 2 ) < 7/(1 - s )> thi s gives p e 1 in L°°(0, T; L 2 7/(7+ 2 ( 1 - a ) 2 )(T 2 )). 
Applying the same procedure as in (|2"T|) again, we obtain 

ll^(Pe _ l)lli<x»(o,T;i 2 7/(7+2(l- s )2)( T 2)) < ||P£||^( 0iT;L V(i- S )(T 2 ))ll^ 7<T ^ £ )ll Loo ( ' T ; L2 ( T2 )) ~^ ' 

Hence, p £ -> 1 in L°°(0, T; W' 1 ' 27 ^ 7+2 ( 1 ~^ 2 )(T 2 )) and in L°°(0, T; L 7 ^ 1_S ) 2 (T 2 )). Repeating 
this argument, we conclude that p e — > 1 in L°°(0,T; L P (T 2 )) for all p < 00. 
For the momentum, we obtain for p > 1 

||p £ U £ — V 0|| ioo ( O)T .£2p/(p+l)(T2)) < |lv / PelU oo (0,T;L 2 P(T 2 ))||A/p7(^£ — V ± 0)||l°o(O,T;L 2 (T 2 )) 

+ 1 1 Per ~~ 1||l oo (O,T;L 2 p/(p+ 1 )(T 2 ))II^" L 0IU oo (O,T;L° o (T 2 )) 
< Cll\/p7(«£ _ V- L 0)||loo( O) T;L2(T 2 )) 
+ C||Pe ~~ 1|U°°(0,T;LP(T 2 ))- 

This shows that p e w e -> V x in L°°(0,T; L q (Y 2 )) for all q < 2. 

Finally, let 7 > 2 and s = 1. Then p £ — > 1 in L°°(0, T; // 1 (T 2 )) and, by the continuous 
embedding H 1 ^ 2 ) ■=— >■ L P (T 2 ) for all p < 00, also p e -> 1 in L°°(0, T; L P (T 2 )) for all p < 00. 
The theorem is proved. 

Appendix A. Local existence of smooth solutions 

The local existence of smooth solutions to the Navier-Stokes-Korteweg system (EQ)-(J5]) 
can be shown similarly as in [21]. We only sketch the proof since it is highly technical and 
does not involve new ideas. First, we rewrite (BJ-Q, setting p = p £ , u = u £ , and e = 1. 
Taking the divergence of flSj) and replacing div<9 t (pw) by (j4j), which has been differentiated 
with respect to time, we obtain 

8l P - + 2^ fi) W f , = - div divO*, « u) - div(^) 

+ 2divdiv(p(p)D(w)) + F[p], 

where F[p] = 2 div(pV(cr'(p) Acr(p))) — 2pa'(p) 2 A 2 p involves only three derivatives. This 
formulation allows one to treat the momentum equation as a nonlinear fourth-order wave 
equation for which existence and regularity results can be applied. In order to derive 
some regularity for the velocity, Li and Marcati [21] assumed that curlw = 0. Then u is 
reconstructed from the problem 

div-u = — (dtp + Vp-u), curlt> = 0, / v(t)dx — u(t). 
P Jt 2 

Theorem 2.1 in [21] gives the existence of a unique solution u e H s+l (T 2 ) to this problem, 
provided that the right-hand side satisfies —{d t p + Vp -u)/p G H S (T 2 ). Actually, Li and 
Marcati replace the right-hand side by —(dtp + Vp-u)/^, where ip solves the mass equation 

d t i/j + ipdivv + u- Vp = 0, t>0, ^(0)=p°. 
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The reason is that this equation can be solved explicitly, yielding strictly positive solutions 
ip. The existence proof is based on an iteration scheme: Given (p p , ip p , u p , v p ), solve 

divup+i = f p (t), curlw p+1 = 0, / v p+1 (t)dx = u(t), 
d t ip p+ i + ip p+1 div v p + u p ■ Vp p = 0, t > 0, ^(0) = P°, 

d ttP P +i ~ ^ A Pp+i + ^pv'H'pf^Pp+i = 9 P {t), t > 0, 

p p+1 (0) = p°, d tPp+1 (0) = -p°div M ° - Vp° • u°, 
d t u p+x + u p+1 = h p (t), 

where f p (t), g p (t), and h p {t) contain the remaining terms (see [211 Section 3] for details). 
The existence of solutions to these linear problems follows from ODE theory and the theory 
of wave equations. The main effort is now to derive uniform estimates in Sobolev spaces 
H k {T 2 ). This is done by multiplying the above equations by suitable test functions and 
assuming that T > is sufficiently small. By compactness, there exists a subsequence of 
(p p ,ipp,u p ,v p ) which converges in a suitable Sobolev space to (p,ij},u,v) as p — > oo. This 
limit allows us also to show that p = if; > and u = v. This shows the existence of local 
smooth solutions under the assumption of irrotational flow curlw = 0. 

Next, we prove the existence of local smooth solutions to the quasi-geostrophic equation 
(IT]). We set p := p(l) > 0. 

Theorem 5 (Local existence for the quasi-geostrophic equation). Let <po G C°°(T 2 ). Then 
there exists T > and a smooth solution <p to (|7|)-(jHJ) for < t < T. 

Proof. The idea of the proof is to apply the theory of linear semigroups. Let p > 2 and 
let A p : W 2,P (T 2 ) — > M, A p (u) = — pAu + u. Then A p is a sectorial operator satisfying 
3?(A) = 1 for all A G o-(A p ), where o~(A p ) denotes the spectrum of A p . Consequently, A p 
possesses the fractional powers A§ for ft > 0, defined on the domain X?* = D(A^). This 
space, endowed with its graph norm, satisfies X 13,1 ' > W k,q (T 2 ) if k — 2/q < 2/3 — 2/p, 
q > p pH Theorem 1.6.1]. Let max{l - 1/p, 1/2 + l/(2p)} < < 1 and set X := X^ p . 
The operator A p generates an analytical semigroup e~ tAp (t > 0) [T71 Theorem 1.3.4], and 
the following estimates hold for alii > [T7| Theorem 1.4.3]: 

||A P e- Mp w|| LP(T2) < Ct-^e- St \\u\\ LPiT 2 h 
\\(e- tAp - i>||z> (T2 ) < Ct p \\ A p v\\ LP{r2) < Ct fi \\v\\ x 

for < S < 1, u G L p (T 2 ), and v EX. 

Next, we reformulate fl7|). Set u = (ft — A</>. Then (J7|) can be written as a system of two 
second-order equations: 

(22) -A(f) + (j) = u inT 2 , i > 0, 

(23) d t u - pAu + u = (V ± cf) ■ V) (0 - u) + p(w - 0) + tt. 



AN ASYMPTOTIC LIMIT OF A NAVIER-STOKES SYSTEM 



17 



We employ a fixed-point argument. Let T > and R > 0. We introduce the spaces 
Y = C°([0,T];X) and B R = {u G Y : \\u-u°\\ Y < R}, where u° = -A0° + 0° G C°°(T 2 ). 
Given u G Y C C°([0, T]; L P (T 2 )), let G L°°(0, T; W/ 2 ' P (T 2 )) be the unique solution to 
( 1221) satisfying the elliptic estimate ||0|| vk 2 -p(t 2 ) < C||w||lp(t 2 )- Then define 

J{u) = e~ tAp u° + / e ( *- s)Ap F(0(s),w(s))ds, where 

JO 

F(0, u) = (V x • V)(0 - u) + fi{u -<f))+u. 

Using the continuous embedding W 2,P {T 2 ) ^ W /1,2p (T 2 ) and the elliptic estimate for 0, we 
infer the estimate 

11-^(07 w)||l°°(0,T;Lp(T 2 )) < ^|M|loo(o,T;W 1 . 2 J'(T 2 )) (l + HwHirc^TjW^Ppr 2 ))) 

< C||w||loo(o,T;X)(1 + \W\\l°°(0,T;X)) 

= C\\u\\y(1 + 

The last inequality follows from the embedding X <^-> H /1,2p (T 2 ) which holds for /3 > 
l/2+l/(2p). 

We show that J maps into .Br and that J : — > Br is a contraction for sufficiently 
small T > 0. Let T > be such that \\{e' tA f - I)u \\ L p^) < CT^\\u°\\ x < R/2. Then, for 
u G B R , 

\\J(u)-u°\\y < sup \\(e- tA - - I)u°\\ LPm 

0<t<T 

+ sup f \\A p e- tAp F(<f)(s),u{s))\\xds 

0<t<T Jo 

sup [\t-s)-^- s ^\\F(4>( S ),u(s))\\ x ds 

^ 0<t<TJo 

if T > is sufficiently small, using that u G Thus J(u) G -Br. In a similar way, we 
show that, for given m,d6 Br, 

|| JN - J(tO||y < -TZTg-dlHk + \Hy)\\u - v\\ Y . 

Again, choosing T > small enough, J becomes a contraction, and the fixed-point theorem 
of Banach provides the existence and uniqueness of a mild solution on [0,T]. 

It remains to prove that the mild solution is smooth. Since /3 > 1 — 1/p, we have 
X W 2 ' p / 2 (T 2 ) and hence u G L°°(0,T; W 2 ^' 2 {T 2 )) C L°°(0,T; W 1,p (T d )). Furthermore, 
V0 G L oo (0,T;iy 1 ' p (T 2 )) C L°°(0, T; L°°(T 2 )) (here, we use p > 2). This shows that 
d t u + Ap(u) G L°°(0,T;L p (T 2 )). Parabolic theory implies that u G L q (0, T; W 2,P (T 2 )) 
for all q < oo. This improves the regularity of to G £ 9 (0, T; W 4 ' P (T 2 )). Hence, 
dtu + A p (u) G L q (0,T; L°°(T 2 )), and a bootstrap procedure finishes the proof. □ 
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